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ABSTRACT 

We report on the discovery of a new Milky Way companion stellar system located at 
(aj2ooo,<5j2ooo) = (22H0 m 43.15 s ,14 o 56'58.8"). The discovery was made using the eighth data re- 
lease of SDSS after applying an automated method to search for overdensities in the Baryon Oscilla- 
tion Spectroscopic Survey footprint. Follow-up observations were performed using CFHT/MegaCam, 
which reveal that this system is comprised of an old stellar population, located at a distance 
of 31.9l^ g kpc, with a half-light radius of — 7.24j^'|g pc and a concentration parameter of 
c = log 10 (r t /r c ) = 1.55. A systematic isochrone ht to its color-magnitude diagram resulted in 
\og(age/yr) = 10.07t°;o3 and [Fe/H] = -l-58t°'°l. These quantities are typical of globular clus- 
ters in the MW halo. The newly found object is of low stellar mass, whose observed excess relative 
to the background is caused by 95 ± 6 stars. The direct integration of its background decontaminated 
luminosity function leads to an absolute magnitude of My = —1.21 ± 0.66. The resulting surface 
brightness is fly = 25.90 mag/arcsec 2 . Its position in the My vs. diagram lies close to AM4 and 
Koposov 1, which are identified as star clusters. The object is most likely a very faint star cluster — 
one of the faintest and lowest mass systems yet identified. 
Subject headings: globular clusters: general — galaxies: dwarf — Local Group 



1. INTRODUCTION 

Recent large surveys such as Sloan Digital Sky Sur- 
vey (SDSS) and the 2-Micron All Sky Survey (2MASS) 
have delivered an enormous amount of data about the 
stellar populations of the Milky Way (MW). These stud- 
ies have probed a new regime in parameter space of 
Milky Way satellites, by significantly expanding the vol- 
ume over which the faintest systems can be detected, 
and have revealed a wealth of new objects. One striking 
discovery from these surveys is the myriad of substruc- 
tures that populate the MW structural components, in- 
cluding stellar streams from disrupting sa tellite galaxies 
and t i dal tails from gl o bular c lusters fe.g iRockosi et all 
' 20021: iMaiewski et al.l 120031: iRocha-Pinto et all 120041 
Newberg et al.l I2010D . It has included the exciting 
discovery of a new class of faint dwarf galaxies (e.g. 
Willman et"all I2005L IBelokurov etall 120061 : flrwin et al.1 
2007tlWalsh etaM om , which may elucidate our under- 
standing of small-scale structure in the dark matter dis- 
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tribution and of galaxy formation at the lowest masses. 
In addition, a handful of very fain t stellar systems have 
been identified in the outer halo (IKoposov et aTl 120071: 
I Belokurov et al.l 120101: iMuhoz et all 120121: iFadelv et all 
120111 ). with somewhat different properties from more 
massive clusters. 

Studies of extragalactic star cluster systems have also 
revealed the existence of diffuse and low surface bright- 
ness stellar systems around luminous galaxies, such as 
Faint Fuzzies (FF) an d Diffuse Star Clu sters (DSC) 
(|Larsen fc Brodid 120001: iPeng et al.l I2006D . Studying 
nearby counterparts of these elusive objects will allow us 
to better constrain their structure, dynamics, and forma- 
tion histories. 

A significant population of MW satellites is very likely 
still left undiscovered (jTollerud et al.l 120081: I Willman I 
120101 ). Future deep surveys, such as the Dark Energy 
SurvejH and the Large Synoptic Survey Telescope El 
will allow very faint systems to be probed much further 
out in the Galactic Halo, and are likely to provide im- 
portant new clues to the formation mechanism of glob- 
ular clusters, to the hierarchical buildup of the Milky 
Way system, to galaxy formation physics in low masses 
systems, and to the abundance and p roperties of the 
small est scale structures in the Universe ()Rossetto et al. I 

unni). 

In this work we present the discovery of a new stel- 
lar system in the MW halo. The object is located 
at (a J20 oo,£j20oo) = (22H0 m 43.15 s , 14°56'58.8"), or 
(l,b) = (75.1735°, -32.6432°), in the Pegasus con- 
stellation. Its standard SDSS object name is SDSS 
J2211+1457, although in this work we choose to call it 
by a shorter name, Balbinot 1. This paper is organized 

10 http://darkenergysurvey.org 

11 http://www.lsst.org 
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as follows. In section 2 we describe the data and meth- 
ods that led to this discovery. In section 3 the follow-up 
observations and data reduction are briefly explained. In 
section 4 we quantify the properties of the object in sec- 
tion 5 and discuss the nature of Balbinot 1 in comparison 
to other satellite systems discovered in the MW. This ob- 
ject is likely a very old stellar cluster, one of the lowest 
mass and lowest surface brightness objects yet detected. 

2. DATA AND METHOD 

The Baryon Oscillation Spectroscopic Survey (BOSS) 
is one of the four surveys that belong to the Sloan Digital 
Sky Survey III (SDSS-III). Its imaging stage is now com- 
plete with photometry in the ugriz system for ~8x 10 7 
sources covering a region of ~ 2000 deg 2 in the south- 
ern Galactic hemisphere was released in January 2011 
as pa rt of the SDSS Data Release 8 (DR8; lAihara et al.l 
2011). So far, no systematic search for MW faint satel- 
lites has been published in the BOSS footprint. Our 
discovery data were taken from the DR8 PhotoPrimary 
view table. A clean sample was defined by selecting those 
objects flagged as reliable stars or galaxies (see detection 
process below). 

To search for stellar systems in a large area of the sky, 
we developed an automated tool c alled FindSat. The al - 
gori thm closely follows the work of lKoposov et al.l ([20081 ) 
and I Walsh et all (|2009f ). There are three main steps in 
the algorithm. First a color-magnitude cut based on stel- 
lar evolutionary models is applied to enhance the pres- 
ence of old metal-poor stellar populations relative to field 
stars. Next, FindSat uses this filtered stellar sample and 
creates a density map on the sky plane. Finally, the den- 
sity map is convolved with a kernel that is the difference 
between two Gaussians, one with the angular size of a 
typical Milky Way (MW) satellite (~ 4 - 8'), and the 
other much wider, used to smooth out any remaining 
large scale structure on the map. 

Source detection is performed on the co nvolved map 
using SExtractor (jBertin fc Arnoutslll996fh To circum- 
vent misidentification due to poor star-galaxy separation, 
we apply identical steps to build density maps of sources 
classified as galaxies by the SDSS reduction pipeline and 
apply SExtractor to them. For a given detected source, 
FindSat then compares the signals in the smoothed stel- 
lar and galaxy density images as determined by SExtrac- 
tor, (S* and Sgai, respectively). Known faint satellites 
tend to fall in the low S ga i and high locus of S ga i 
vs. S* plane. We use this region to select our new satel- 
lite candidates. We also compare the position of these 
new overdensities with those of known objects found in 
Abell, Globular Clusters, and NGC catalogs, to remove 
any known object that may share the same locus of the 
Sgai vs. S* plane. 

For the BOSS region we ran FindSat using color- 
ma gnitude cuts based o n four Padova evolutionary mod- 
els |Girardr~et~irj [2002) with ages of 8 and 14 Gyr and 
Z = [0.001,0.006]. We tuned our CMD filter to search 
search for objects at five distance moduli (m — M) = 
[18.0, 19.0, 20.0, 21.0, 22.0]. Visual inspection of the de- 
tected candidates was carried out, revealing Balbinot 1 as 
a faint overdensity of blue sources consistent with them 
being stars. However, deeper and higher spatial resolu- 
tion images were required to investigate the nature of 
Balbinot 1. 
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Figure 1. Bottom panel: the isodensity contours built using the 
list of sources detected by Terapix. The contour levels are for 
[8,16,24,32,40] stars/aremm 2 . The dashed box shows a 1' X 1' 
region for which we show the combined g band image in the top 
panel. 



3. FOLLOW-UP OBSERVATIONS 

In the first semester of 2012, we obtained follow-up 
observations of Balbinot 1 using MegaCam, which is in- 
stalled in the Canada- France-Hawaii- Telescope (CFHT). 
MegaCam is a mosaic 36 CCDs, each with 2048 x 
4612pix. The total field of view (FOV) is 0.96 x OMdeg, 
with a pixel scale of ~ 0.19 " Jpix, 

We designed the observations to reach the magnitude 
of the main sequence turn-off (MSTO) of a 10 Gyr stellar 
population at 150 kpc distance with S/N ~ 10. Expo- 
sure times of 2800 s in the g band and 3800 s in the 
r band were required. The total exposure time was di- 
vided into 6 dithering exposures for each passband to 
avoid scattered light from bright stars, saturation, and 
also to cover the smaller gaps between the CCDs. The 
observations were carried out under photometric condi- 
tions, seeing was always below 0.8", and the airmass be- 
low 1.15. 

The basic reduction of the images (overscan and bias 
subtraction, and flatfielding) was done by the CFHT 
team. The images were geometrically corrected, regis- 
tered and co-added by Terapix, which uses the Astro- 
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matic toolkitPl. Therefore, source detection is based on 
SExtractor. Point spread fun ction (PSF) photometry 
was carried out using DAOPHOT (|Stetsonlll9M ). The input 
coordinates list for the stars measured using DADPHDT was 
the one provided by Terapix. This list presented less spu- 
rious detections than the one from DAOFIND. Photometric 
calibration was achieved using bright and non-saturated 
(19 < gsDSS < 21) stars in the observed field, whose 
instrumental magnitudes were then compared to those 
from SDSS. We found 536 matches which we used to 
find a calibration equation composed of a photometric 
zero point and a color term. The equation coefficients 
were found by means of a least square fit using sigma 
clipping rejection. Our mean calibration residuals are of 
0.024 in the g band and 0.019 in the r band. 

The observed MegaCam field shows lit tle differential 
redden ing. Using the dust maps from ISchlegel et al.l 
(1998) we find a value of E(B-V) = 0.060 at position 
of Balbinot 1 and a standard deviation of 0.0014 across 
the whole field. 

Thus, differential reddening is not likely to contribute 
significantly to our analysis. Nonetheless we choose to 
correct for reddening using the value of E(B-V) at the 
position of each star in the Mega Cam field. We adopt 
Rv = 3.1 and the coefficients from lCardelli et all (|1989[ ) 
to compute the extinction in the CFHT passbands. 

Figure [T] shows a stacked g band image centered on 
Balbinot 1 (top panel). We also show isodensity con- 
tours of our detected sources (bottom panel). A clear 
overdensity of stars is seen on both panels of the figure 
very close to where the FindSat candidate was originally 
identified. 

4. RESULTS 

From Figure[T]it is evident that there is a concentration 
of stellar-like objects in the observed region. To confirm 
its stellar nature we must analyze its CMD. In Figure 
[2] we show the extinction corrected (g — r) x g CMDs. 
In the top left panel we show the CMD based on the 
SDSS discovery data. On the top right panel we show 
the CMD from our follow up CFHT images. In both 
cases the stars are restricted to a radius of 150" from 
the visual center of Balbinot 1. To discard Balbinot 1 as 
a density fluctuation of field stars we also show (bottom 
left panel) the CMD for a ring well away from the object 
center covering the same area on the sky. By comparing 
these panels we not only clearly confirm the excess of 
stars around the position of Balbinot 1 but also conclude 
that the distribution of these stars on the CMD plane 
is consistent with a simple stellar population (SSP). An 
obvious MSTO is seen at (g - r) ~ 0.21 and g ~ 21.4, 
in connection to an RGB stretching up to (g — r) ~ 0.6 
and g ~ 18.5. 

The bottom right panel in Figure [2] shows the MSTO 
region in detail. We also sho w the best fitting Padova 
isochrone (jGirardi et al.ll200"2h . It corresponds to an age 
of 11.7ig g Gyr and an abundance of Z = 0.0005, or 

[Fe/H] = — 1.58±g:?| at a distance of d© = 31.9±i £ kpc. 
In the same panel we show the effect of varying the dis- 
tance modulus within its fit uncertainties. The isochrone 
fit was performed by means of a x 2 minimization of the 

12 http://www.astromatic.net/ 
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Figure 2. Top left panel: (g — r) X g CMD for SDSS measured 
sources inside a radius of 150" from the visual center of Balbinot 
1. Top right panel: the CMD for the same region as the previous 
panel but showing sources from our CFHT photometry. The white 
triangles show the object ridge line. Bottom left panel: CMD for 
sources located in a ring far away from the object center. The 
area of the ring is the same as the area of the inner region. In the 
bottom right panel we show zoomed in version of the top right panel 
around the MSTO region; we also show the photometric errors for 
each source. The best fit isochrone is shown as a solid line. The 
isochrones with upper and lower limits of distance modulus are 
shown by the dashed lines. 



CMD distance between the cluster ridge line and the 
isochrone set. Our model grid goes from 9.8 to 10.12 
in log(agefyr) with a step of 0.01, and from 0.0001 to 
0.004 in Z, with steps of 0.0001. We adopt a step of 
0.01 in distance modulus and a range from 16.0 to 18.0. 
The best model is the one that minimizes x 2 - The un- 
certainties are derived from models that have values of 

Y 2 = Y 2 ■ +1. 

A. Amm 1 

We stress that the best fit isochrone and its associated 
parameters are somewhat dependent on the stellar evo- 
lution model and on the model grid. In fact, our quoted 
metallicity uncertainty is comparable to the Padova grid 
resolution for this parameter. Also, the discrepancies 
among different evolutionary models are much larger 
than our quoted uncertainties. For a comprehensive ap- 
proach to this problem we refer to iKerber fc Santiago! 
(2009]). 



4.1. Radial profile and half-mass radius 

To quantify the size and concentration of Balbinot 1 
we built the radial density profile (RDP) by counting 
stars in annuli around its center. To increase contrast 
relative to the background, only stars within 2cr in color 
away from the best fit isochrone were selected, where 
a is the mean photometric error at a given magnitude. 
The data was also cut at g < 24 to avoid photometric 
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Figure 3. The radial density profile for Balbinot 1 with la error 
bars in the y direction. The error bars in the x direction are the 
bin sizes. The solid line shows the best fit King profile. The dashed 
vertical line shows the position of the core radius. The dot-dashed 
line shows the half-mass radius, and the dotted line shows the 
corresponding plummer profile. 

incompleteness. The result is shown in Figure [3] where a 
very peaked distribution of stars is visible. 

In order to avoid the binning complications that arise 
when trying to fit a RDP we chose a Maximum Likeli- 
hood (ML) approach. The method we u s e her e follows 
closely the one described in iMartin et al.l (|2008f ). 

The ML fit was carried out using stars in the range 
0.0 < R < 12' and fitted both a King and a Plum- 
mer model. This limit was chosen to avoid the large 
CCD gap of MegaCam. The best fitted models are 
shown as the solid and dotted lines in Figure EJ respec- 
tively. For the King profile we derive r core = 0.22to'og' 
and r -tidal — 7.85+ 3 ; 40 '. And the Plummer profile 

fit yields r s — 0.60^°'^' and a total number of stars 
M* = 94.7±5.9. The uncertainty in each parameter is es- 
timated using the likelihood-ratio approximation, where 
likelihood values near the maximum approaches a chi- 
squared distribution. The uncertainties are taken as the 
90% confidence level of the chi-square distribution. For 
the Plummer profile we also kept the center of Balbinot 1 
as a free parameter. This allowed us to find the best cen- 
ter and its uncertainty. The best center was then used in 
the King profile fit. Notice that the upper bound of the 
tidal radius uncertainty is outside the range of R avail- 
able in our observations. For simplicity we have adopted 
the upper bound as the largest radii available, that is 12'. 

For the Plummer profile the half-mass radius is eas- 
ily obtained from the relation 77, = 1.305r s , yielding 
the value of r h = 0.78 tg;^'. For the King model we 
made an estimate of the half-mass radius as follows. We 
first integrated the profile from zero to the limiting ra- 
dius and subtracted the expected number of background 
stars. The result of this operation yields the total num- 



ber of observed stars that should belong to Balbinot 1 
N oos . We then compute the half- mass radius as the ra- 
dius which contains N i> s /2 stars, again taking care to 
subtract off the expected background sources. We obtain 
r/j = 0.69^Q 25'- The two estimates of agree within the 
uncertainties. Since the Plummer model is widely used in 
most studies of satellites, we choose to use the half-mass 
radius estimated from it. 

4.2. Total luminosity 

The total luminosity of Balbinot 1 is dominated by 
bright red giant branch (RGB) stars and possibly a few 
red clump (RC) stars, although Figure [2] does not reveal 
any obvious RC star in the inner 150" region. The small 
number of evolved stars in the RGB and, most especially, 
in the RC, indicates that Balbinot 1 is a low mass stellar 
system, more consistent with being a star cluster than 
a dwarf galaxy. Given its low stellar mass, to properly 
estimate the total luminosity of Balbinot 1 one must take 
into account the background (and foreground) star con- 
tamination. 

We first built the observed luminosity function (LF) of 
Balbinot 1 by counting stars as a function of magnitude 
inside a circle of 150" radius from the cluster center and 
in a ring with inner border at 300" and outer at 550". 
We subtracted the area weighted star counts in the ring 
from the area weighted counts in the inner circle. This 
results in an observed LF for Balbinot 1 already decon- 
taminated from field stars. The adopted ring is 10 times 
larger than the cluster circle to minimize fluctuations in 
the subtracted field contamination. To build the LFs 
we again applied a magnitude cut of g < 24.0 to avoid 
incompleteness . 

The total magnitude is obtained by direct integration 
of the background decontaminated LF. Using the uncer- 
tainty on each bin of the LF and a bootstrap method 
we derive the uncertainty on the total magnitude. The 
final value is My — —1.21 ± 0.66 for Balbinot 1. Cou- 
pling the My value with the r% estimate yields a surface 
brightness of fly — 25.90 mag/arcsec 2 . The value of My 
of Balbinot 1 is comparable to the value obtained for 
Muhoz 1. This luminosity is also comparable to tha t of 
SEGUE 3 (|Belokurov et alj|2010b fFadely et al.ll2011| ). 

5. SUMMARY AND DISCUSSION 

In this paper we report on the discovery of a new stellar 
system in the MW halo, found the SDSS-III/BOSS foot- 
print in the Southern Galactic hemisphere. Its confirma- 
tion as a genuine stellar system required deep follow-up 
imaging from CFHT. By means of a theoretical isochrone 
fit, we derived a heliocentric distance of 31. 9*t{. 6 kpc , 
an old age of 11. 7^8 Gyr, anc ^ a metallicity of [Fe/H] 
= —1.58^0 i3- We also found that a King profile provides 
a good description of its structure; the best fit profile has 
a core radius of r c = 2.04+ ' 56 pc 1 , a limiting radius of 
r t = 72.84^31^5 pc, and a projected half-mass radius of 
r/j = 7.24j^ 29 P c obtained using a Plummer profile. We 
carefully estimate the object total luminosity by means of 
direct integration of the background decontaminated LF. 
With the aid of a statistical bootstrapping we estimate 
the uncertainties on the total absolute magnitude of Bal- 
binot 1 leading to the final value of My — —1.21 ± 0.66. 
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Table 1 

Summary of the derived parameters for Balbinot 1. 



Parameter 


Value Unity 


"J2000 
<5j2000 
d Q 

My 

log(agefyr) 
[Fe/H] 


22 : 10 : 43.15 ±0.3 h:m:s 
14 : 56 : 58.8 ±2.0 ° : ' : " 
31.9±i-g kpc 
-1.21 ±0.66 mag 
10.07too3 dex 
— 1.58±g;?| dex 


Plummcr 


r s 

S f>9 

M* 
rh 


0.6+JJ;" arcmin 
1.38to'o5 stars/arcmin 2 
94.7 ±5.9 
0-78±qH arcmin 


King 


r c 
rt 


0.22+JJ;^ arcmin 
7.85t3 4 o arcmin 
1.32+^;^ stars/arcmin 2 
162.76 ±5.0 stars/arcmin 2 
°- 69 -0.25 arcmin 



In Table 1 we present a summary of the parameter de- 
rived for Balbinot 1. 

The total number of stars and absolute magnitude of 
Balbinot 1 suggest that it is a star cluster. Its size is 
larger than that of most clusters i n the Galactic system , 
either open or globu l ar, or in M31 (ISchilbach et al.l f2006: 
Ivan den Bergbll20l"ol 120111 ). However, it falls close to the 
median radius when compared to outer halo clusters in 
the Galaxy. Its size is more typic al of the diffuse star 
clusters found bv lPeng et al.l ()2006[ ) in early- type galax- 
ies, but again with several orders of magnitude difference 
in terms of luminosity. It is an extremely low luminosity 
cluster; In fact, Balbinot 1 seems to be one of the faintest 
and lowest surface brightness old stellar systems found so 
far in the MW. Only Mu noz 1 has a lower surface bright- 
ness (|Munoz et al.ll2012l ). and only Munoz 1 and Segue 3 
have lower absolute magnitudes. Its location in the lumi- 
nosity vs. size diagram places Balbinot 1 close to other 
systems identified as low lumin osity outer halo clu sters, 
including Koposov 1 and AM4 (jFadelv et al.H2011l) . 

The stellar clo ud that lies closest to Balbinot 1 is the 
Hercules- Aquila (jBelokurov et al.ll2007h . The cloud is lo- 
cated at I ~ 50° although its extension is poorly known, 
specially in the southern galactic hemisphere. However 
the heliocentric distance of the cloud is 10 — 20 kpc, mak- 
ing it unlikely that Balbinot 1 is associated with this halo 
structure. 

The small size of the object, and the lack of any clear 
evidence for complex stell ar populations make it un likely 
that it is a dwarf galaxy (|Wilhnan fc Strader|[2012l ). We 
note that this locus is at the confluence of the branches 
filled by classical gl obular clusters and M W dwarfs, as 
shown for instance in iMcConnachid (|2012l ). Spectra from 
individual stars are being obtained and will allow mea- 
surement of metallicity spread and line-of-sight veloci- 
ties, which may help determine the dynamical mass of 
Balbinot 1 and constrain its stellar population. 

1 The scale transformation from sky to physical values is 9.28 
pc/arcmin 
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